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The influence of gravitational lensing on the spectra of 
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ABSTRACT 

We consider the influence of (miUi/micro)lensing on the spectra of lensed QSOs. We 
propose a method for the observational detection of microlensing in the spectra of 
lensed QSOs and apply it to the spectra of the three lensed QSOs (PG 1115+080, 
QSO 1413+117 and QSO 0957+561) observed with Hubble Space Telescope (HST). 
We find that the flux ratio between images Al and A2 of PG 1115+080 is wavelength- 
dependent and shows differential magnification between the emission lines and the 
continuum. We interpret this magnification as arising from millilensing. We also find 
that the temporal variations in the continuum of image C of QSO 1413+117 may be 
caused by microlensing, while the temporal variation observed in QSO 0957+561 was 
probably an intrinsic one. 
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1 INTRODUCTION 

Gravitational macro and microlensing are well known phe- 
nomena that have been widely discussed in the litera- 
ture (e.g., Schneider, Ehlers & Falco 1992; Zakharov 1997; 
Narayan & Betelmann 1999; Wambsganss 2001, Claeskens 
& Surdej 2002 and references therein). The influence of mi- 
crolensing on the spectra (continuum and spectral lines) of 
lensed QSOs has been investigated mainly theoretically (Ne- 
miroff 1988; Schneider & Wambsganss 1990; Wambsganss & 
Paczynski 1991; Popovic et al. 2001a; Abajas 2002; Lewis & 
Ibata 2004). Some of the observational effects have also been 
presented in several papers (e.g., Lewis et al. 1998; Mediav- 
illa et al. 1998; Wisotzki et al. 2003; Wucknitz et al. 2003; 
Chelouche 2003; Metcalf et al. 2004; Richards et al. 2004). 

The influence of microlensing on the line profile of 
lensed QSOs was initially investigated by Nemiroff (1988) 
and Schneider & Wambsganss (1990). Recently Wandel et 
al. (1999) and Kaspi et al. (2000) using reverberation tech- 
nique discovered that the broad line region (BLR) size was 
smaller than that assumed in the standard AGN model. Fol- 
lowing this discovery, the influence of microlensing on the 
line profile of lensed QSOs in the UV/optical wavelength 
band, was revised and discussed in several papers (Popovic 
2001a; Abajas 2002; Lewis & Ibata 2004). Similar investi- 
gations have been performed in the X-ray spectral region 
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(Popovic et al. 2001b; Chartas et al. 2002,2004; Popovic et 
al. 2003; Dai 2003). 

Abajas et al. (2002) noted that the influence of mi- 
crolensing on the shapes of spectral lines can be very strong 
for high ionization lines arising from the inner part of the 
broad emission line region. The magnification depends on 
the size of the BLR and we therefore expect that distortions 
of spectral lines will not be observable in all microlensing 
events. Using the BLR size estimation for NGC 5548 (Pe- 
terson & Wandel 1999) and the relation between the ra- 
dius of the BLR and optical luminosity, Rblr ~ 
A = 5100 A, given by Kaspi et al. (2000), Abajas et al. 
(2002) identified a group of 10 gravitational lensed systems 
in which microlensing of the BLR could be observed. 

Observations and modeling of microlensing of the BLR 
are promising, because the study of the variations of the 
broad emission line shapes in a microlensed QSO image 
could constrain the size of the BLR and the continuum re- 
gion (e.g., Lewis & Ibata 2004). Observations of microlensing 
of the profiles of Fe Ka spectral lines in three multi-imaged 
QSOs were recently presented in Chartas et al. (2002,2004) 
and Dai et al. (2003). The observed variations of the Fe 
Ka lines are in good agreement with theoretical predictions 
(Popovic et al. 2001a; Popovic et al. 2003). Observations of 
the effect of microlensing on the profile of broad UV emis- 
sion lines (C IV and S IV/OIV) were recently reported by 
Richards et al. (2004). 

The detection of microlensing from the distortion of 
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the spectra of lensed QSOs is more complicated than in- 
ferring its presence from photometric observations. First of 
all, larger telescopes arc needed to obtain the high signal- 
to-noise (S/N) ratio spectra, preferably with 2D spectropho- 
tometry that can provide simultaneous spectra of all images. 
In addition, the intrinsic variability of lensed QSOs needs 
to be taken into account when comparing spectra of differ- 
ent images. Gravitational microlensing can also distort the 
spectrum of the continuum. Specifically, it is expected that 
microlensing is wavelength dependent (e.g., Wambsganss & 
Paczynski 1991; Lewis et al. 1988; Wisotzki et al. 2003) since 
the size of a continuum emission region of the accretion disk 
depends on the wavelength band. 

According to the standard model of AGNs, a QSO con- 
sists of a black hole surrounded by a (X-ray and optical) 
continuum emitting region probably with an accretion disk 
geometry, a broad line region and a larger region that can 
be resolved in several nearby AGN that usually is referred 
to as the narrow line region (e.g., Krolik 1999). 

Variability studies of QSOs indicate that the size of the 
X-ray emission region is of order lo'^*"^® cm (e.g., Char- 
tas et al. 2001; Oshima et al. 2001, Dai et al 2003). An Fe 
Kq fluorescence line detected in several AGN near 6.4 keV 
is thought to originate from within a few 10s of gravita- 
tional radii (e.g., Tanaka et al. 1995; Fabian et al. 1995). 
This line is thought to be a fluorescence line of Fe due to 
emission from a cold or ionized accretion disk that is illu- 
minated from a source of hard X-rays originating near the 
central object. Aspects of this emission can provide a probe 
of strong gravity near a black hole. 

Gravitational lensing is in general achromatic (the de- 
flection angle of a light ray does not depend on its wave- 
length), however, the wavelength dependent geometry of the 
different emission regions may result in chromatic effects. 
Studies aimed at determining the influence of microlensing 
on the spectra of lensed QSOs need to account for the com- 
plex structure of the QSO central emitting region. Since the 
sizes of the emitting regions are wavelength dependent, mi- 
crolensing by stars in the lens galaxy will lead to a wave- 
length dependent magniflcation. The geometries of the line 
and the contirmum emission regions arc in general different 
and there may be a variety of geometries depending on the 
type of AGN (i.e., spherical, disk-like, cylindrical, etc.). 

On the other hand, in some cases the potential of the 
lens galaxy may be perturbed by small satellite galaxies 
or globular clusters (hereafter millilcnsing) . These pertur- 
bations will add additional complexity to the magnification 
function. 

To investigate the effects of microlensing, a method is 
needed that can be applied to the observed spectra. A tech- 
nique that is often used (in cases where flux losses from slits 
are significant) relies on measuring the equivalent widths of 
spectral lines (Lewis et al. 1998; Wisotzki et al. 2003). One 
limitation of this approach is that it does not account for 
microlensing of the continuum component that may signifi- 
cantly affect the line equivalent widths. 

The aim of this paper is to discuss the influence of gravi- 
tational (micro/milli)lensiug on the spectra (continuum and 
spectral lines) of lensed QSOs. First in §2 we outline a model 
that includes the geometry of the line and continuum emis- 
sion regions, as well as the magnification due to the lens 
galaxy and the stars in the lens galaxy. Based on this model 



we propose a method for detecting microlensing events in 
the lines as well as in the continuum components of QSOs. 
In §3 we applied the method to a sample of lensed QSOs 
observed with HST, and in §4 we present our conclusions. 



2 THE INFLUENCE OF LENSING ON QSO 
SPECTRA 

We assume that the emitting region of a QSO is geomet- 
rically complex, and that the geometries and sizes of the 
continuum and the line emitting regions are different, i.e., 
that the unlensed surface brightness, /(A; X, Y) (in the case 
of X-ray I(\;X,Y) I{E;X,Y), where E is the emitted 
energy), can be written as a function of wavelength as, 

7(A; X, Y) = r (A; X, Y) + /^(A; X, Y), (1) 

where I''{\;X,Y) and l'"{X;X,Y) are the surface bright- 
nesses of the continuum and the lines, respectively and X,Y 
are the coordinates of the emission region in the source 
plane. The magniflcation of the images depends in part on 
the projected mass distribution of the lens galaxy. We as- 
sume that the macro- magnification (A(X,Y)) is complex and 
can be represented as function of the X,Y coordinates in the 
source plane. Under these assumptions an observer will de- 
tect the following magnified intensity of an image (i), 

9i(A;X,y) =9J(A;X,y)-h9f(A;X,y) 

where 

S>i'^(A;X,y) =7'=-^(A;X,y) •^(X,y) (2) 

In general one can expect an energy dependent magni- 
fication of the images of a lensed QSO resulting in different 

magnifications of the line and the continuum emission re- 
gions (as it was observed in the case of Q2237+0305, see 
Mediavilla et al. 1998). 

Taking into account that the sizes of the continuum and 
the line emitting regions may be a function of wavelength, 
we write the fluxes as, 

*?'''W = J r'^(A;X,y)-A(X,y)dS, (3) 

where S is the projection of the QSO emitting region in 
the source plane onto a plane perpendicular to the line-of- 
sight of the observer and dE = dXdY. As one can see from 
Eq. (3), in general, the fluxes, $i(A) and $c(A), depend on 
the geometry of the source as well as on the magniflcation 
caused by the lens galaxy. 

If we assume there are no millilensing and the sizes of 
the emission regions are signiflcantly smaller than the size 
of the macro-caustic, then we expect the magnification of 
the contiimum and the line fluxes to be the same, i.e., 

$=(A) = F=(A) • A., $f(A)=F^(A).^,. (4) 

where Ai is the macro-magnifications of the image i and 
^"^'■^(A) = J^r'^(\,X,Y)dT.. 

The inference of a microlensing event (MLE) from the 
analysis of the spectral lines and the continuum spectral 
components of lensed QSOs can be complicated by several 
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effects that may include intrinsic variability of their fluxes 
and their spectra and the fact that intrinsic variability be- 
tween different images of a lensed QSO will be time-delayed. 

A possible solution to this problem is to obtain ob- 
servations of all the images of a lonscd QSO at different 
epochs, preferably with the time intervals between observa- 
tions coinciding with the time-delays of the system. With 
such observations one can compare the flux of each com- 
ponent as a function of wavelength for different epochs, as 
well as the flux ratios between different images observed at 
a given epoch. We define the observed intensities of fluxes 
A and B as ^a{X) and $s(A) at time to and ^'aW and 
$'b(A) at time ti, and the time-delay between the images as 
At^t-i- to. 

In principle, images of a lensed quasar are always being 
microlensed by stars in the foreground galaxy. This results 
in the continuous magnification as well as demagnification of 
the lensed quasar images. In the quiescent periods, when the 
magnification is low, the magnification map is smooth and 
one should not expect any detectable differences between the 
spectral energy distributions of images. In the case where 
magnification caustics are in the vicinity of the source the 
differential magnification of the continuum is expected to 
be significantly stronger and detectable (e.g., Wambsgauss 
& Paczynski 1991). Consequently, hereafter we consider only 
high-magnification (Amle) microlensing events. 

If a microlensing event (with a strong amplification) is 
present in only one of the images (e.g., A at time ti), then 
the fiux of this image will vary with time, and the magnified 
flux due to microlensing can be written as: 



$UA) 



Ai / [r(A;X,y) + 7^(A;X,y)]- AMLEiX,Y)dTi (5) 



where E is the projected surface of the macrolensed image 
at the source distances and Amle is the amplification due 
to microlensing. The flux ratio {Ra',a = ^a(A)/$a(A)) of 
the same image at two different epochs separated by the 
time-delay is: 



Ra'a 



(A; X, Y) + J^(A; X, Y)]Amle{X, Y)dT, 
[F=(A)+F^(A)1 



(6) 

and the flux ratio between different images observed at the 
same epoch is 



Ra 



/j,[/=(A; X, Y) + 7^ (A; X, Y)]Amle{X, y)dS 



A2 [F=(A) + F^(A)] 

We define as an indicator of variability (cr): 



<t{\)^ Ra',a- Rb' 



(7) 



(8) 



In general a = Q indicates no variability. We note the rare 
case where tr = and variability is present. This condition 
can occur when the fiuxes of both images A and B change by 
the same factor between epochs to and ti. The probability 
of such an event is negligible. If cr ^ 0, then we investigate 
the source of the variability. 



If microlensing of only one image (e.g. A) is present 
and there is no other source of variability, then the average 
magnification caused by microlensing is. 



Amle{^) 



[r {X;X,Y)+I^{X;X, Y)]Amle {X, Y )dE 

From Eqs. (5-7) we write the average magnification 
caused by microlensing as 



KA) = 



'&a(A) _ A2 '^''aW 
$a(A) Ai *^(A)' 



(9) 



where Ai and A2 are the macro-magniflcations of images A 
and B, respectively. 

Eq. (9) is only valid for the case where microlensing is 
occurring in one component. For the case where intrinsic 
variability (cr/) is present we have 



ffi — Ra\a — -t-Ra'b' 7^ 0, 



(10) 



Also, ai will be different from zero if both images are 

microlensed at the same time. This may be checked by com- 
paring the spectra of two images observed at two different 
epochs separated by the time-delay. If the relation 

_^A' _ Ai 

holds we conclude that intrinsic variation is present, other- 
wise if it is not satisfied we conclude that both images are 
microlensed. 

On the other hand, taking into account that the total 

flux is the sum of line and continuum fluxes (see Eq. (1)), 
the average magnification of only the line is 

-jL $A'(A)-$^,(A) _ 

AmleK'''^) — 



^a(X) - *i(A) 

A2 $A'(A)-$^,(A) 
^1 *B,(A)-*^,(A)' 



(11) 



and in the case of millilensing one can expect that mag- 
nifications in the line and in the continuum are different 
{A^ ^ A°) as well as that these magnifications are wave- 
length dependent. 

In Table 1 we point out several cases where our method 
can be used. There arc more complex cases not listed in 
Table 1 such cis the combination of milli and microlensing, 
microlensing and intrinsic variability and microlensing of 
both images. Observations in several epochs separated by 
the time-delay are useful in distinguishing between these 
complex cases. 



3 OBSERVATIONS: SEVERAL EXAMPLES, 
RESULTS AND DISCUSSION 

Our method of analyzing the spectra of lensed quasars out- 
lined in §2 can be applied accurately when simultaneous 
observations of spectra of images of a lensed QSO at two 
different epochs separated by the time-delay are available. 
In many cases the time-delay of a system is unknown. In 
such a situation our method can still be applied to estimate 
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Table 1. A list of cases where our proposed method can be used to infer the presence of millilensing, microlensing and 
intrinsic variability. 



The case Description 



I) <T = and Rab = R-a'b' = const. 


No variability or millilensing present 


Ila) (7 = 0, Rab = Ra'b' = fW and/or jt A" 
lib) (7 = 0, Rab = Ra'b' = const, and A^ ^ 


There is no variability, but millilensing or differential extinction 
might be present (sec e.g., Wucknitz ct al. 2003), if Rab = /(A) 
cannot be explained by extinction — > millilensing is present 


III) (7 7^ 0, (7j 7^ 0, and Ra'b = ^ 


Intrinsic variability is present 


IVa) a^O,ai = 0, Ra,b ¥= fW and Ra',b' = fW 


Microlensing of one image at time ti 


IVb) (7 ^ 0, (7/ = 0, Ra,b = /i{A), Ra',b' = /2(A) 
and Ra'a = 1 (or Rb'b = 1) 


Microlensing of one image at times to and ti 
(timescale of microlensing is longer than time delay) 



the origin of the magnification. We note that the applica- 
tion of our method is limited by tlie S/N of the component 
spectra, e.g., in the case of low S/N, non zero values for the 
variability indicator sigma and chromatic effects may only 
be detected in cases where there are large flux magnifica- 
tions caused by microlensing or other effects listed in Table 
1. Therefore high S/N spectra are needed to apply our model 
accurately. 

Here we apply our method to several observations of 
lensed quasars performed with HST. By inspection of the 
HST database we found only three lensed quasars with si- 
multaneous observations of all images at different epochs: 
QSOs Q0957-f561, PG 1115-f080 and Q1413-I-117. In this 
section we compare the spectra between different images as 
well as between the same image at different epochs using 
the method described above. The STARLINK DIPSO soft- 
ware package (Howarth et al. 2003) was used to analyze the 
spectra. 

3.1 Q0957+561 

Images A and B of Q0957-I-561 have been observed with HST 
at several epochs. In the HST archive we found observations 
of both images made with the FOS G270 and obtained in 
1995 (Jan 26, Oct 20, Nov 04, Nov 17, and Dec 14) that 
we used for our analysis. We also found observations of only 
image A of Q0957-F561 with STIS performed on 1999 April 
04, and observations of only image B on 2000 June 2, and 
only image A on 2000 September 8, but these observations 
have not been used in our present analysis. A detailed de- 
scription of all observations and the analysis of the spectra 
can be found in Michalitsianos et al. (1997), Dolan et al. 
(2000) and Hutchings (2003). 

A first inspection shows that images A and B have sim- 
ilar spectra (practically the same), with the exception of the 
spectra observed on 1995 January 26 (see Figure 1). We ap- 
plied the method described above to all observed spectra of 
Q0957-I-561 obtained in 1995, comparing the spectra of im- 
ages A and B from different epochs. From observations taken 
within the period October-December we did not detect any 
significant difference between the spectra that would indi- 
cate a microlensing event or intrinsic variability. A slight but 
not statistically significant difference exists in the variability 
indicator a (see Figures 2a, 3a and 5a). A significant differ- 



ence between images was only detected in spectra observed 
on 1995 Jaimary 26, (Figures 2b, 3b and 5b). In Figure 2 we 
present the variability indicator between observations per- 
formed on October 20 and December 14, 1995 (Fig. 2a) as an 
example. The variability indicator a (similar to the case of 
all spectra observed within the period October-December), 
docs not significantly differ from zero and taking into ac- 
count the noise in the spectra we conclude that significant 
variability is not present in this time interval. On the other 
hand a comparison between the spectra observed on 1995 
January 26 and spectra taken at other epochs shows signif- 
icant differences. As an example the variability indicator a 
between observations performed on Jan 26 and Dec 14, 1995 
is shown in Figure 2b. Variability is clearly present in the 
eontirmum as well as in the Lya line (the location of the line 
is indicated with a dashed line). To investigate the nature 
of this variability we determined the fiux ratio using Eq. 6. 
As shown in Figure 3, the flux ratio for the observed spectra 
in the period October-December is practically one, whereas 
the flux ratio of the spectra observed on 1995 January 26, is 
^1.8 times higher than the ratio observed during the period 
October-December. 

We do not flnd any signiflcant wavelength-dependence 
of the variability indicator a, the flux ratio Ra,ai between 
the same image at different epochs or the flux ratio Ra,b 
between different images for the same epoch. 

Moreover, in Figure 4 we show that at three different 
epochs within the October-December period the flux ra- 
tio Ra,b is consistent with the magnification resulting from 
macrolensing. Using the optical magnitudes of images A and 
B of Q0957-I-561 available from the Castles database^ we ob- 
tained a flux ratio of Fa ~ l.OS^s . A comparison of the flux 
ratios of Figure 4 also shows that on 1995 January 26, the 
flux of the continuum of image B was signiflcantly larger 
than that of image A. 

We use Eq. (10)^ to investigate the nature of the vari- 
ability in Q0957-f561. As shown in Figure 5, the variability 
(cTj 7^ 0) seen on January 26, 1995 is most likely intrinsic. 

Here we briefly discuss the absorption line detected 



^ http: / / cfa-www.harvard.edu/castles/ 

^ Here we approximate the amplification ratio as = ^B. ^ 
Rb,a 
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Figure 1. The spectra of Q0957+561 observed with HST for 
three different epochs. The soUd Une and dots indicate the 
spectra of images A and B, respectively. The flux is given in 
erg cm~^sec~^A~^. 
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Figure 2. The variability indicator a = (I'^/Ia) — {Ig/Is) cal- 
culated for the spectra of Q0957+561 for the epochs; a) 1995 
December 14 and 1995 October 20 (case I); b) 1995 January 26 
and 1995 December 14 (case II). The vertical dashed line indicates 
the location of the Lya line. 
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Figure 3. The flux ratio between the spectra of image A of 
Q0957+561 for the epochs: a) 1995 December 14 and 1995 Octo- 
ber 20 (case I); b) 1995 January 26 and 1995 December 14 (case 
II). 
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Figure 4. The flux ratio between images A and B of Q0957-f561 
observed for the listed epochs. 



blue-wards of Lya originating from a damped Lyman al- 
pha (DLA) system at a redshift of ~ 1.391 (Michalitsianos 
1997). A comparison of the flux ratios between images at all 
epochs indicates that the absorption line is stronger in im- 
age A than in imago B (sec Figure 4). This is in agreement 
with the previous analysis reported by Michalitsianos et al. 
(1997). Such differences in the absorption between images 
indicates inhomogeucity in the DLA system that may arise 
from the different paths of the rays of the two images (see 
e.g. Chelouche 2003). 

To summarize, our analysis of the spectra of Q0957-I-561 
taken in January, October and December of 1995 indicates 
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Figure 5. crj, an indicator of intrinsic variability of Q0957+561 
spectra is shown for the epochs: a) October and December 1995 
and b) January and December 1995. 



that: i) the strong magnification detected in the continuum 

and in the Lya lino is most likely caused by intrinsic vari- 
ability of the QSO, ii) during this variation image A has a 
stronger line flux than image B and image B has a stronger 
continuum flux than image A, iii) the flux ratio of the con- 
tinuum is not wavelength-dependent. 



3.2 PG 1115-1-080 

PG 1115+080 is a lonscd QSO with a redshift of 1.72. The 
time-delay between images B and C of 23 ± 3.4 days was 
initially measured by Schochter et al. (1997) whereas a later 
analysis by Barkana (1997) resulted in a value of 25±2 days. 
For the purpose of applying our method we ignore the time- 
delay between images Al and A2. Impcy ct al. (1998) had 
reported an anomalous flux ratio between images A2 and 
Al. Their observations ruled-out differential extinction and 
microlensing as explanations for the anomalous flux ratio, 
because of a lack of variability. They conclude that since 
an expected flux ratio between images A2 and Al of ~ 1 
is a generic feature of the large-scale potential near a fold 
caustic (from a simple lens model), only a potential per- 
turbation intermediate between that produced by stars in 
the lens galaxy (microlensing) and by the overall galaxy can 
explain the anomalous flux ratio, and that the potential of 
PG 11154-080 must be perturbed either by a satellite galaxy 
or a globular cluster, i.e. that millilensing is present. 

Our search of the HST database found an observation 
of PG 1115+080 made on 1995 June 07, which covers the 
wavelength range from 850 A to 1200 A (rest wavelength) 
and observations made on 1996 Januaxy 21 (image Al) and 
on 1996 January 24 (image A2) that cover the wavelength 
range from 850 to 1750 A (rest wavelength). These observa- 
tions were performed with the G270H and G400H gratings. 
More details describing these observations can be found in 
Bechtold et al. (2002). 

First, we calculated the indicator of variability (a) and 
intristic variability (aj) between observations made on 1996 
January 21 (image Al) and on 1996 January 24 (image A2) 
for the wavelength range 850 A - 1200 A. As shown in Figure 
6a, the variability indicator is not significant, but it appears 



that a low level of variability may be present (cr ~ 0.1). In 
Figure 6b we show that weak intrinsic variability may be 
present (a/ « 0.05). However, both indicators of variability 
are negligible with respect to the level of the rms noise (~ 
0.3). 

To further investigate the significance of variability, we 
determined the flux ratios between the same images at two 
epochs. As shown in Figure 7, variability is not significant. 
Specifically the continuum fiux of image Al increased by 
about 3%, and that of image A2 decreased by about 7% be- 
tween July 1995 and January 1996. (see Figure 7), however, 
both variations are insignificant relative to the level of the 
rms noise. 

We have also calculated the flux ratios between images 
Al and A2 for the July 1995 and January 1996 observations 
(see Figures 8). As shown in Figure 8 the flux ratios tend to 
be slightly different for different epochs, not only in magni- 
tude, but also in the shape of the wavelength dependence. 
This result should be taken with caution because of the high 
level of noise, especially in the spectra obtained on 1996 Jan- 
uary. We also find that the emission lines in the spectrum of 
image A2 were stronger than those in image Al during the 
July 1995 and January 1996 observations. This difference in 
lino strength is more obvious in the fiux ratio between im- 
ages Al and A2 of PG 1115+080 for the observation made 
on 1996 January when plotted in the spectral range of 850 
A - 1750 A(see Figure 9). The N IIIA990.98, O IVA1033.82, 
Lya, Si IVA1396.76 and C IVA1549 lines are stronger in im- 
age A2 than in imago Al. Wc note that the flux ratio as 
defined in Eqs. (6) - (7) includes both line and continuum 
components. Consequently, the difference of the intensity of 
the emission lines in images Al and A2 of PG 1115+080 as 
shown in Figures 8 and 9 may arise in part from the different 
magnifications of the underlying continua in images Al and 
A2. 

We assumed that the relative macro-magnification of 
images Al and A2 is achromatic (i.e., case II, Table 1). To 
search for possible differences between the magnification of 
the spectral lines and continuum, we scaled the spectra ob- 
served in January 1996 of both images on the same con- 
tinuum level. In Figure 10a we compare the flux vs. wave- 
length of images Al and A2 of PG 1115+080 after scal- 
ing the continuum of image A2 to that of Al (multiply- 
ing the flux of the image A2 continuum by the function 
R{A1,A2) = 2.75 - 5.9 ■ lO""* ■ A). We notice that the spec- 
trum of image A2 has stronger lines (with respect to the 
local continuum) than those of image Al. This is expected 
if we have microlensing or millilensing of the continuum. But 
this difference cannot be explained by microlensing since the 
indicator of variability is small with respect to rms noise 
(cr ~ 0.1). The difference between the continuum and line 
amplifications is more likely the result of millilensing. Here 
we should mention that long-term microlensing may produce 
a similar effect (case IVb), even with the indicator of vari- 
ability being negligible in this case, we cannot completely 
exclude long-term microlensing. 

If we assume that the indicator of variability is negli- 
gible (a ~ 0) and millilensing is present (case (2) of §2.1), 
then we expect millilensing to cause different magnifications 
of the continuum and the emission line components, as indi- 
cated in Eq. (3) . Since the line and continuum emission com- 
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ponents arc additive (Eqs. 1-3), wc subtracted the contin- 
uum in both images (Eq. 11) and compared the line fluxes. 

A comparison of the line fluxes in the spectra of images 
Al and A2 with subtracted continua is shown in Figure 10b. 
A small difFerences in these fluxes can be due to uncertain- 
ties in the subtracted continua, but wo can conclude that 
the A1/A2 flux ratio of the lines is {A1/A2)l ~ 1. This ra- 
tio corresponds to the expected value of ~ 1 if no micro or 
millilensing were present because the images are symmetri- 
cally arranged near a fold caustic (see Schneider et al. 1992) . 
On the other hand, the ratio of A1/A2 continuum flux is 
wavelength dependent and in the observed wavelength range 
(from 850 A to 1750 A) can be approximated by the function 
R{A1, A2) continuum « 2.75— 5.9- 10"'' -A, i.e. the amplification 
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Figure 9. The flux ratio of images Al and A2 of PG 1115-1-080 
in the wavelength band of 850 - 1750 A 



in the continuum is changing from (j41/A2)continuuin « 2.3 
at 850 A to (Al/A2)continuum » 1-75 at 1750 A. 

Since millilensing is similar to microlensing, we may ex- 
pect the line profiles of images Al and A2 to be different 
(Popovic et al. 2001, Abajas et al. 2002). In Figure 11 we 
compare the Lya and C IVA 1549 line profiles between im- 
ages Al and A2. We also show a comparison between the 
Lya and C IVA 1549 line profiles from the spectrum of im- 
age A2 since the line emitting regions of different emission 
lines may be different (e.g., Wandel et al. 2000; Kaspi et al. 
2000). 

For the purpose of these comparisons, we normalized 
the fluxes of the Lya and CIVA1549 lines to one and con- 
verted wavelengths to velocities : A — > X = (A — Ao)/Ao. As 
shown in Figure 11 there are no significant differences be- 
tween the line profiles of the same lines in the spectra of im- 
ages Al and A2. Those comparisons suggest that millilensing 
does not affect the emission line region of PG 1115-1-080. 

In Figure 11c, we note that the red wing of the Lya 
line is stronger than the C IV red wing. We propose that 
this difference is likely caused by the contribution of the 
NVA1240 fine. We note that such a strong NVA1240 doublet 
(Fjv V ^ O.SFLy) is unusual for QSOs, with a more typical 
value expected to lie in the range Fn v ~ (0.1 — 0.2)Flj, 
(e.g., Laor 1994). 
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Figure 10. The flux (given in erg cm^^sec^^A^^ units) of image 
Al of PG 1115-1-080 (dashed line) compared to that of image A2 
of PG 1115-1-080 (solid line): a) after scaling the spectra of both 
images to the same continuum (the continuum level of image Al), 
and b) after subtraction of the continuum in both images. 



Our investigation of the spectra of images Al and A2 
of PG 1115-1-080, confirms the anomalous spectral distribu- 
tion and relative magnifications of these two images, first 
noted by Impcy ot al. (1998). The main results of our anal- 
ysis of the spectra of PG1115-I-080 arc as follows: i) We 
only detect very weak variability in images Al and A2 of 
PG 1115-1-080 between two observations; ii) the differential 
magnification of the line and the continuum components as 
well as wavclcngth-dcpciidcnt amplification may be caused 
by millilensing. This is consisted with an earlier proposal by 
Impey et al. (1998) that suggested the presence of significant 
perturbations to the main lens galaxy of the PG 1115+080 
lens system by a satellite galaxy or globular cluster, iii) 
Millilensing affects the continua of Al and A2 in observed 
wavelength range {F{Al)c/F{A2)c « 1.75 - 2.3, expected 
~ 1), but it does not appear to significamtly affect the line 
profiles of the images {F{A1)l/F{A2)l « 1). 



3.3 QSO 1413+117 - Cloverleaf 

The spectra of the images of QSO 1413+117 were observed 
in the spectral ranges of 900-1350 A and 1300-1800 A and 
an analysis of these spectra was presented in Monier et al. 
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Figure 11. A comparison between : a) the Lya line profiles of 
images Al and A2 of PG 1115+080; b) the C IV line profiles of 
images Al and A2 of PG 1115+080; c) the Lya and C IV hne 
profiles of image A2 of PG 1115+080 (the vertical dashed line 
shows the position of the N vA1240 line). 



(1998). They found that the emission spectral lines show an 
intrinsic absorption component and that the emission lines 
in the spectrum of image D arc significantly weaker than in 
the other three images. 

For our analysis we used the observations performed on 
1993 June 23, (images A and D) and 1993 June 27, (images 
B and C) covering the spectral range 1300 - 1800 A and 
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Figure 12. The flux ratio (Ri.j) of different images of QSO 1413+117 observed on the same date. Left panels: Ri j versus wavelength 
ranging from 900 A to 1350 A. Observations were made on 1994 December 24. Right panels: Rij versus wavelength ranging from 1300 
A to 1800 A. Observations of images A and D were performed on 1993 June 23 and of images B and C on 1993 June 27. Prom top to 
bottom, we show the flux ratios between the following images: A and B, A and C, B and C, A and D, B and D; C and D. The solid and 
dashed lines indicate the best fits of Ri^j = a + b ■ \ and i?, j = const, respectively. 
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December 24, 1994 (all images) covering the spectral range 
900 - 1350A. The number and separation in time of these 
observations are not ideal for an accurate application of our 
method, however, as we will show we can still use our method 
to discuss the presence of a Icnsing event. 

In Figure 12 we show the flux ratios as a function 
of wavelength (for both spectral ranges) of all images of 
QSO 1413+117. We fit the flux ratios with models of the 
form R = a + b ■ X and R = const. 

An inspection of the flux ratios indicates that: 

i) during the first epoch (1994 December 24) only the ra- 
tios R{A,C) and R{B,C) have a significant non-zero slope, 
i.e., they arc wavelength-dependent. This might also be the 
case for R{C,D), but due to large residuals in the emission 
lines and the high level of noise of the ratio this result should 
be taken with caution. In the other cases the ratios tend to 
be constant. 

ii) the flux ratios R{A,B), R{A,D), and R{B,D) show 
no significant differences between all epochs observed, but 
flux ratios containing image C are different between epochs, 
i.e., R{C, D), R{B, C), and R{A, C) show significant changes 
of the fitted slope between epochs. 

iii) the spectra of image D for all observations show 
weaker UV emission lines than those of the remaining images 
as was previous noted in Monier et al. (1998). 

From these observations we conclude that significant 
variability only of image C is present, and considering that 
the ratios R{A,C) and R{B,C) tend to be wavelength- 
dependent, we suspect that microlensing of image C was 
likely present during the observations. 



4 CONCLUSIONS 

In this paper we investigated the influence of gravitational 
lensing on the spectra of lenscd QSOs. Starting from the 
assumption that the magnification due to a lens is in general 
a complex function, (i.e., the presence of globular clusters 
or satellite dwarf galaxies in the lens galaxy may introduce 
perturbations to the potential of the principal lens galaxy as 
noted by Impey et al. 1998) and that the line and continuum 
emitting regions are different in sizes and in geometries, we 
found that the magnification of the spectra of the different 
images may be chromatic (as was noted in Wambsganss & 
Paczyiiski 1991; Lewis et al. 1998; Wisotzki et al. 2003 and 
Wucknitz ct al. 2003). 

Here we briefly summarize several conclusions that arise 
from Eqs. (1-11) in §2: 

i) in general the complex geometries of the omitting re- 
gions of QSOs and the complex potentials of the lens galax- 
ies need to bo taken into account when determining their 
influence on the observed spectra of the continuum and line 
emission; 

ii) microlensing events load to wavelength-dependent 
magniflcations of the continuum that can be used as indica- 
tors of their presence; 

iii) the average magnification due to microlensing may 
be determined from the flux ratio of an image observed at 
two different epochs separated by the time-delay and from 
the flux ratio between different images observed at the same 
time; 

Here we propose a method to infer the presence of mi- 



crolensing, millilensing and intrinsic variability in one or 
both of the images of a lensed quasar. The method is out- 
lined in the following steps: 

(a) Obtain simultaneous observations of the spectra of 
the images of a lensed quasar (at least) at two different 
epochs separated by the time-delay. 

(b) Calculate the variability parameters a and ai using 
equations (8) and (10), and determine the flux ratios Ra'a, 
Rb'b, Rab, Ra'b', and Ra'b- 

(c) Use Table 1 to determine the possible presence of mi- 
crolensing, millilensing or intrinsic variability of the quasar. 

Using our proposed method we investigated the influ- 
ence of microlensing on the FOS spectra of three lensed 
QSOs. We presented an example of the chromatic effect in 
our analysis of the FOS spectra of PG 1115-1-080. In partic- 
ular, the application of our lensing method to this system 
indicates that the line and continuum components of the 
lensed QSO spectra of PG1115-I-080 have different magnifi- 
cations during the HST observations. 

The application of our method to investigate the effect 
of lensing on the spectra of lensed QSOs may be useful not 
only for constraining the unresolved structure of the central 
regions of QSOs, but also for providing insight to the com- 
plex structure of the lens galaxy and possible substructure 
in the lens galaxy. 
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